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Abstract 

Studies of magnetic induction in von Karman swirling flows have so far linked the time-averaged 
induced magnetic field to the structure of the mean flow. They have evidenced the Omega and 
Parker mechanism generated respectively by the fiow differential rotation and helicity, which un- 
derly the Duddley and James l| dynamos. Using an array of Hall probes we study here the 
dynamical regime. In the experimental flow, turbulence is fully developed and large fiuctuations 
are observed in the magnetic induction processes. We find that the large scale turbulent fiuctu- 
ations have different characteristics when induction results from the differential rotation or from 
the dynamics of the stagnation point in the mid plane of the von Karman fiow. Symmetry con- 
siderations indicate that the dynamical fiow spends half of its time away from the time-averaged 
structure. The consequences of these observations for dynamo experiments are discussed. 

PACS numbers: 47.65. +a.47.27.-i.91.25.Cw 
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I. INTRODUCTION 



The von Karman flows generated in the gap between counter-rotating impellers has been 
considered by several groups as a possible candidate in the search for a laboratory 

demonstration of a homogeneous fluid dynamo that would be less constrained than the Riga 
and Karlsruhe designs 0,01. In these experiments, liquid sodium is often used because its 
electrical conductivity is high (half that of copper, at 120°C), while its density remains low, 
of the order of that for water. Its magnetic diffusivity, as for all metals, is many orders 
of magnitude larger than its hydrodynamic viscosity. The flow Reynolds number needs 
to be very large in order for non-linearities to develop in the magnetic induction. As a 
consequence, the hydrodynamic flow is very turbulent, and many questions arise concerning 
the influence of turbulence on the bifurcation threshold and the dynamics in an eventual 
saturated regime. 

This problem is very complex, and many studies so far have focused on the dynamo 
capacity of the average flow engineered in the von Karman (VK) geometry. Here, 'average' 
has the meaning of 'time-average'. A time-averaged flow fleld (U)(r) is derived from mea- 
surements in water prototypes (for convenience and because velocimetry methods are scarce 
in opaque, high temperature flows of liquid metals 0,0) as: 

1 '■^ 



(U)(r) = - lJ{r,t)dt, (1) 

where T is a time much longer than the one characteristic of the forcing of the flow (esti- 
mated, for instance, as the period of rotation of the driving impellers). This stationary flow 
profile - no longer a solution of the Navier-Stokes equation - is then usually inserted into a 
kinematic numerical solver, in which the induction equation 

9tB = Vx ((U) xB) + AAB , (2) 



(A is the magnetic diffusivity) is solved with the velocity field kept constant in time j9|,|l0|, 



These studies have shown the possibility of dynamo action in the average VK fiows, and the 
underlying induction processes have been measured in sodium and Gallium fiows ll^ and 
analyzed in details ^,14, 3|- It showed that the helicity and differential rotation present 



in the von Karman mean flow cooperate to generate a self-sustained dynamo. 

The existence of a kinematic dynamo threshold for the VK flows, together with the possi- 
bility to bring its value within experimental reach (in terms of power requirements) has mo- 



tivated the sodium experiments in Maryland, Cadarache and Wisconsin ( Q] and references 
therein). However, it has also been recognized early that VK flows at high Reynolds 
numbers have strong fluctuations. As a result the instantaneous mean flow structure can 
differ signiflcantly from the time-averaged flow. This has motivated investigations regarding 
the influence of noise on the dynamo bifurcation. For instance, it has been proposed that 
the VK flow may lead to an intermittent dynamo burst of dynamo activity would 

occur, triggered by transient flow structures that are most efficient at generating magnetic 
induction. 

Induction measurements in the presence of an externally applied field have been made 
in the VKS experiment jj]. They have demonstrated the existence of fluctuations in the 
magnetic induction due to small scale turbulence and also generated from the instationarity 
of the large scales. For example, it has been observed that the local fluctuations of induction 
measured during a time interval r increase logarithmically with r j^. This flnding can be 
linked to the existence of a 1// scaling in the low frequency domain of the magnetic induction 
spectrum, a feature that has been also observed in numerical models j^l and in the dynamo 
regime of the Karlsruhe experiment [2^. 

The purpose of this article is to explore further the above ideas. We relate the fluctuations 
in time of magnetic induction to large scale fluctuations in the geometry of VK flows. We 
use a VK flow of liquid Gallium, and we study the magnetic induction in the presence of 
an externally applied fleld. Besides being relatively easy to handle in the laboratory, the 
magnetic Reynolds number in Gallium flows reaches values of order 1: induction effects are 
measurable and the non-linearities are not yet strong. The amplitude of the induced fleld is 
often an order of magnitude smaller than the amplitude of the applied fleld, and, in the quasi- 
static limit 3|; the induction equation Q is reduced to V x ((U) x Bo) + AAB = 0, with Bq 
the applied fleld. When it is uniform in space, the induced fleld solves AAB = — (Bq ■ V)(U) 
so that it is directly related to the velocity gradients in the flow (although the solution is not 
local; the overall topology of the velocity gradients and the boundary conditions do enter in 
the solution of this Poisson equation). 

The induction measurements are made using a new probe made of a line of magnetic 
fleld sensors. It samples simultaneously one component of the magnetic fleld at several 
location within the flow. The paper is organized as follows: in the next section, we describe 
the experimental set-up and we characterize our multi-probe magnetic measurements. The 
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results are detailed in section III and their implications regarding dynamo action in VK 
flows are discussed in section IV. 

II. SET-UP AND MEASUREMENTS 

A. Flow. 

Our experiments are carried out in the setup sketched in FigO The flow is produced by 
the rotation of one or two disks inside a stainless steel cylindrical vessel flUed with liquid 
Gallium. The cylinder radius R is 97 mm and its length is 323 mm. The disks have a 
diameter equal to 165 mm and are fltted a set of 8 blades with height 10 mm. They are 
separated by a distance H = 203 mm. The disks are driven by two 1 1 kW AG-motors which 
provide a constant rotation rate in the interval Q G [0.5, 25] Hz with a stability of about 
0.1%. The system is cooled by a set of coils located behind the driving disks; the experiments 
are made with the flow kept at a temperature in an interval between 42°G and 48°G. The 
fluid is liquid Gallium (density p = 6.09 x 10^ kg m~^) whose electrical conductivity is 
a = 3.68 X 10^ ohm~^ m~^. Its kinematic viscosity is u = 3.1 x 10~^ m^ s~^. The integral 
kinematic and magnetic Reynolds numbers of the flow are deflned as Re = 2t[B?VL/v and 
Rin = 27ifj,QaR'^fl. Values of R^ up to 5 are achieved, with corresponding Re in excess of 
10^ 

B. Induction measurements. 

Two sets of coils (see Fig. Q]) are used to apply an external magnetic fleld Bq, either 
parallel or perpendicular to the axis of the cylinder. The locations of the coils is such that 
the configuration is close, but not strictly equal, to a Helmholtz geometry. Variations of 
applied field intensity over the fiow volume are of the order of 18% of the mean, for axial 
or transverse applied fields. Its magnitude Bq is less than 100 G, so that the interaction 
parameter, 

N = aBlL/pU = aB^2R/p27rRn ~ lO'^ , 

is quite small. One can safely neglect the back-reaction of the Lorentz forces on the velocity 
field. 
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FIG. 1: Set-up : (a) schematics of the von Karman mean (time-averaged) flow geometry for counter 
rotating disks at equal rate 0,. The rotation is opposite in each half of the cylinder, leading to 
a strong differential rotation and the recirculation loops in the flow create a stagnation point: 
the poloidal flow converges in the mid-plane and is drawn towards the disks on either side, (b) 
experimental arrangement: coils and magnetic measurement probe. 

Magnetic induction measurements are performed using a multi-sensor probe. Eight Hall 
sensors form a linear array which can be inserted inside the flow to yield simultaneous 
measurements of the magnetic field along a line, i.e. measurements of magnetic induction 
profiles. In the experiment, the probe array is inserted radially into the flow, in the 
mid-plane between the driving disks. The magnetic fleld is sampled at 8 locations between 
1.5 and 8.5 cm from the rotation axis. The spacing between the Hall sensors is equal to 
1 cm, equal to the overall diameter of the probe - small compared to magnetic diffusive 
scale in the Gallium flow {^m ~ 2R/RJ^^'^ > 6 cm). According to the orientation of the 
probe, one obtains a radial proflle of one component of the magnetic fleld, Bx{ri,t) or 
Bz{ri,t) (for i = 1...8) — the axis are shown in Fig. ^ The Hall sensors are single axis 
elements from Sentron (ISA-IM), with a sensitivity of 0.03 V/G and a frequency range 
from DC to 10 kHz. We have calibrated these probes by comparison with a temperature 
compensated Hall probe connected to a Bell FW-9953 gaussmeter. Their temperature 
dependence is small, of the order of 0.01 G/K in the temperature range of interest i.e. 
between 40 and 50°C. In typical measurements, the signals from the 8 sensors are recorded 
for durations between 150 and 480 seconds, using a National Instrument PXI-4472 digitizer 
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at a rate of 1000 Hz with a 23 bits resolution. 




f(Hz) 

FIG. 2: Spectrum obtained with one Hah sensor, compared to a measurement made using a 
calibrated Bell gaussmeter. Flow driven by the rotation of a single disk, at 17 = 10 Hz. The main 
peak is at 8.5 Hz because the flow rotation rate is lower that that of the driving disk. 



We compare in Fig. |21 a typical frequency spectrum computed using a time series from a 
single sensor, compared to the equivalent measurement made with the Bell gaussmeter. In 
this case, the flow is set into motion by the rotation of only one of the driving disk, as in 
ref. this choice is made in order to have a steady flow with little large scale fluctuations. 
The two curves are in excellent agreement. The dynamical range spans 80 dB. Our in- house 
magnetic line array has thus characteristics that are equivalent to those of the commercial 
single probe. 



III. RESULTS 



A. Average induction profiles 



We present here the analysis of mean induction effects which will be useful for our dis- 
cussion of the fluctuations - the reader is also referred to |l2|. 

We label Bq the applied field and B(r, t) the field induced by the flow motion. We write 
the total magnetic field as Bq + (B) + b, where (B) is the time average of the induced field, 
and b the fluctuating part ((b) = 0). By definition, the k-th component of time-averaged 
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induced field at location Tj is 



{B,{rj)) = ^ [ dtBk{r,,t). (3) 




f (fr) 

FIG. 3: (Top) Time-averaged induction profile {BQ){rj) for an axial applied field Bq^z- (left): 
measurements for Rm = 1, 2, 3, 4, top to bottom (the symbols are the measurements and the solid 
line is the result of the numerical simulation based on the time-averaged flow), (right): schematics 
of the Omega effect responsible for the induction. 

(Bottom): corresponding figures for {Bz){rj) measurements in the case of a transverse applied field 
Bq^x- In this case the induction is due to the stretching and compression of the field lines by the 
axial flow. 

In Fig. 3(a,b) we show the induction profiles {Bg){rj) and {Bz){rj) measured when an 
external field Bq = Bqz is applied, directed along the axis of the cylinder. As detailed 
in previous studies (e.g. Q]), the azimuthal component (Bg) is mainly generated by the 
counter rotation induced by the driving impellers on each side of the mid-plane, via the 
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Omega effect [2 



X{A{B))e = -Bo^A{Ue) • (4) 



On the other hand, when one probes the induced field along the axis, the main effect is 
due to the axial stretching cause by the motion of the ffuid towards each driving disks (the 
stretching component of the stagnation point in the mid plane). 

XA{B,) = -Bq,MUz) ■ (5) 

These mechanisms are sketched in figure 3. We stress that these profiles are obtained from 
simultaneous sampling in space of the magnetic field, and not from a series of independant 
measurements in which a single probe is successively located at a set of radial positions, as 
previously reported [J, [l^ . 

We then compare the measurements to a numerical solution of the stationary induction 
equation: 

= V X ((U) X (Bo + (B))) + AA(B) . (6) 



In the numerica 



measurements 19 1 and the induction equation is solved using the iterative method described 

Q.J 
(the scheme is based on finite difference computations, and thus allows for 

the implementation of realistic electromagnetic boundary conditions). Because Rm remains 

small (less than 5), we have found that in each case it was sufficient to compute the induction 

up to order 3 in Rm. In Fig. 3, the profiles obtained numerically are shown as solid lines. 

One observes an excellent agreement with the experimental measurements. It indicates that 

the mean induction is very well accounted for by the time-averaged flow — equation ©. 

Using all combinations for the directions of the applied and measured induced magnetic 

fields, we have observed that the induction is dominated by two main mechanisms: the 

differential rotation generated by the counter rotating disks (the toroidal part of the flow) 

and the existence of a stagnation point in the mid plane (the poloidal component). For an 

applied field parallel to the rotation z-axis (case in Fig. 3), the toroidal velocity gradients 

induce a toroidal magnetic field (the Omega effect) and the poloidal velocity gradients 

generate a stretching of the applied field lines. When the field is applied transverse to the 

rotation axis (e.g. along the x-axis) the poloidal flow produce a compression of the applied 

field lines while the toroidal flow generates a ^-component associated to the connection 



computation, the time-averaged velocity field is taken from experimental 
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of the magnetic field lines induced along the y-axis on each side of the mid plane. This 
effect is quite sensitive to the electromagnetic boundary conditions; we have called it the 
BC-eflFect Q- 

Regarding the mean (time-averaged) value of the magnetic induction, our main finding is 
that the average induction coincides with the induction predicted from the average flow — 
at least for the magnetic Reynolds numbers reached in this Gallium experiment. One has 

(B) = induction from (U) , (7) 

i.e., the mean induced field effectively solves V x ((U) x Bq)) + AA(B) = 0. We will 
report elsewhere a detailed analysis of the induction effects due to the turbulent small scale 
fluctuations Q|; its main finding is to confirm the above statement: the small scale turbulent 
fluctuations do not contribute to the mean induction. 

B. Fluctuations 

1. General considerations 

As pointed out in the introduction, due to the very small value of the Prandtl number, 
the flow is fully turbulent. The Reynolds number is larger than 10^, even for the relatively 
moderate values of the magnetic Reynolds number reached here. The local velocity fluctua- 
tions reach 35% of the mean speed, and so do the induced fields, as shown in figure 4(a,b). 
The fluctuations are large and occur over a very broad range of time scales. The fastest 
fluctuations are associated with the shearing of the applied field by the turbulent small scale 
motion Komogorov's scaling for the inertial range of motions together with the use of 
Taylor's hypothesis predict an /"^^^^ for frequencies larger than that, fi, of the forcing. This 
is due to the fact that in our experiment Rm is of order one, so that the magnetic dissipation 
scale iu is of the order of the flow size. Our measurements - figure 4(c) - show a steeper 
slope, of the order of —4.6 for driving disks fitted with blades. The Kolmogorov prediction 
— 11/3 was observed in the flow generated by flat rugose disks [2^. The slope of the spectra 
in the high-frequency (dissipative) region does not depend upon which component is being 
probed, nor on the rotation rate of the disks. 

Another observation is that the magnetic induction also fluctuates in a broad range of 
long time scales, as indicated by the behavior of the power spectra for frequencies lower than 
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the disks rotation rate Q - Fig. 4(c). In addition, it is clear in figure 4(c) that the lower 
end of the frequency spectra does depend on which induction mechanism is at work. In the 
case of an axial applied field, the fluctuations in the toroidal component of the induction are 
related to fluctuations in the differential rotation. The low frequency part of the spectrum 
is close to an /~°'^ behavior. The fluctuations of the induced field in the direction of a 
transverse applied field are linked to the dynamics of the stagnation point in the mid plane. 
In this case the low frequency part of the spectrum has stronger fluctuations (the spectrum 
is close to a /^^'^ behavior). The fluctuations of induction are associated to fluctuations in 
the position of the stagnation point. Concerning this slow dynamics, we have also varied 
the disks rotation rate between 5 Hz and 20 Hz without detecting any noticeable change in 
the above features; the spectra as in figure 4(c) collapse when rescaled by the rotation rate 
and the fields rms amplitude. 

■1" , „ 
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FIG. 4: Fluctuations for local measurements. Counter-rotation of the disks at = 10 Hz. (a) 
Bx induced component for an applied axial field i?o,z) (b) induced component for an applied 
transverse field -Bo,xi (c) Corresponding time spectra (black line : signal in (b); blue line : signal in 
(a)). The curves have been shifted vertically for clarity; (d) Power spectral density of the pressure 
recorded at the flow wall in the mid-plane. 

The slow dynamics in the induction traces back to the evolution of the hydrodynamic 
flow. This is evidenced by measurements of pressure fluctuations at the wall; one observes a 
1// behavior in the same range of frequencies where the magnetic fields shows a long-time 
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dynamics - figure 4(d). However, the magnetic measurement is quite sensitive because is 
probes the velocity gradients selectively. The poloidal and toroidal components are inde- 
pendantly probed by chosing the direction of the applied field and the particular component 
of the induced field under inspection. 

2. Profile fluctuations 

We analyse below the fluctuations of induction 'profiles' which we define as the set of 
measured magnetic field values {Bk{ri,t), ...Bk{r8,t)} with {ri,...,rs) radial positions at 
distances between 1.5 and 8.5 cm from the axis of rotation, and the applied field along 
direction i. The choice of the {i, k) couple determines which component (toroidal or poloidal) 
of the flow is probed. We stress again that we do not address here the question of the 
contribution of the small scales of turbulence. The recorded signals Bk{rn,t) are low-pass 
filtered with a corner frequency equal to 3 times the disks rotation rate Q. 
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FIG. 5: Examples of fluctuations of profiles. Counter rotation of the disks at = 10 Hz; applied 
field along the z-axis, Sq.z = 24 G. (left) transverse induced field due to differential rotation; (right) 
axial induced field due to the stretching effect. 

Examples of the evolution in time of such profiles are shown in figure 5. The profiles 
are shown at time intervals equal to the period of rotation of the disks, and one observes 
significant fluctuations with respect to the time-averaged curves. In addition, the points 
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in the profiles tend to vary as a block. This is evidenced by computing the correlation 
function, between, say the first element in the array, and the progressively more distant 
ones. The resulting variation is shown in figure 6(a), where one observes that the maxima 
of the correlation functions decrease very slowly from 0.9 for the second element down to 
about 0.4 for the last one. These values do not change as the disks rotation rates are varied 
from 5 Hz to 20 Hz. For an estimation of a typical correlation length, one notes that the 
curve can be fitted by an exponential function exp — (r/ro), with tq ~ R, the radius of the 
cylindrical vessel. This confirms that the fluctuations in the induction are a global feature 
rather than a local one. In addition, we have observed that the correlation functions are 
symmetric and peaked around the zero time-lag value. This behavior is also evidenced in 
figure 6(b) which shows a space-time diagram of the profile dynamics: fluctuations are felt 
simultaneously at all points, i.e. with a response time of the order of 




FIG. 6: (a) Evolution with the distance to the first element in the array, of the maximum (reached 
at r = 0) of the correlation function Cij(r) = {B{ri,t)B{rj,t + r)). Counter rotation at 17 = 
5, 10, 15, 20 Hz, corresponding to Rm values between 1 and 4. (b) Space time diagram, at = 15 Hz. 
Measurements for counter rotating disks, with Bq^x = 48 G, induced field along the same direction 



3. Distance to the mean 



Let us compute a global distance between the mean flow induction profile and a realisa- 
tion. Using the £2 norm, we define 



N=8 



i=l 



(8) 
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We find that the time averaged value of E are found to be of the same order of magni- 
tude as the average induced field, that is {E(t)) ~ (B). In units of the applied magnetic 
field, we measure (E) ~ 0.07 Bq for situations in which the maximum of the induced field 
(max,r^(i?)(rj)) is also of the order of 0.1 Bq. These observations indicated that instantaneous 
induction profiles differ significantly from the time averaged computation. In addition, we 
also observe that the distance E has large fluctuations away from its mean value. This is 
evidenced in figure 7(a) which shows the probability density functions V{E) for induction 
measurements probing the toroidal and poloidal flows. The curves are wide, with events 
that span several standard deviations. The fluctuations are found to be larger in the case of 
the induction due to differential rotation than when due to the stretching by the stagnation 
point in the mid-plane. 




FIG. 7: Probability density function of the fluctuations of the C2 distance E{t) to the mean 
induction profile. Counter rotation at = 10 Hz, corresponding to Rm = 2. 



As in our experiment the intrinsic magnetic Reynolds numbers (defined by the ratio 
max{B)/BQ 26[) reached in our experiment is always less than 1, we believe that these 
fluctuations mirror the variations in the flow velocity gradients. However the geometry of 
these changes is lost with the use of a global £2 norm. In the next section we analyse further 
the profile fiuctuations, with the purpose of trying to quantify how the induction deviates 
from the one expected from the mean flow. 
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4- Polynomial analysis 



Because of the strong correlation in the signal measured by successive elements in the 
magnetic array, the induction profiles are smooth and are well described by polynomials 
of order three. We do not claim here that the fitting functions are the actual solutions 
to the induction equation. We use the polynomials as a way to take advantage of the 
symmetries associated to the von Karman flow. For instance, the mean velocity is reflection 
- symmetric about the mean plane, and the azimuthal velocity is zero on the rotation axis. 
As a result, the induction due to differential rotation in the presence of an axial applied field 
(the Q effect) has a zero average in the mid plane at r = 0. Deviations from this value can 
then be associated with a symmetry breaking in the mean flow pattern. We thus write an 
instantaneous profile as 

{Bk{r,t)} = ao{t) + ai{t)r + a2{ty + asity (9) 

and we study in the sequel the evolution of the coefficients aj{t). 

• Mean flow polynomials 

We first present the coefficients (aj) for the mean profiles, and their evolution with the 
disks rotation rate. From the results presented in section 3.1 we expect that they are well 
accounted for by the structure of the mean flow, i.e. (aj) = aj(from (U)). 

Let us first consider the evolution of in the case of an axial applied field (along 
z). As detailed in section 3.1, induction in the mid plane is dominated by the twisting 
of the magnetic field lines by differential rotation. Thus, (oq) should be null because of 
axisymmetry and (oi) should increase linearly with the disk (counter) rotation rate. This 
is indeed observed in figure 8. We have no simple interpretation for the coefficients a2 and 
as which mainly ensure that the induced field vanishes at the outer cylinder (because of the 
insulating boundary condition). 
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FIG. 8: Polynomial coefficients for the mean induction profile. Axial applied field Bq^z = 24 G. 
The measured profiles and corresponding polynomial fit are for Rm = 1,2,3,4. (top): induction 
due to differential rotation; (bottom): axial stretching of the applied field. 
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FIG. 9: Polynomial coefficients for the mean induction proffie. Transverse applied field, i?o,x = 
48 G. The measured proffies and corresponding polynomial fit are for Rm = 1,2,3,4. (top): 
induction due differential rotation; (bottom): transverse compression of the applied field in the 
mid-plane. 
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Similarly, when ones probes the ^-induced field the main effect is due to the stretching of 
the applied field lines. Since the applied field is uniform, one expects a non zero component 
(ao), which increases with the disk (counter) rotation rate - as confirmed in figure 8(c). The 
other contribution is in the 02 coefficient, since here the axisymmetry requires that the odd 
terms be null. 

The situation is equivalent when the applied field is transverse Bq = Bqx. The average 
field induced along x in the mid plane comes from the compression of the applied field by 
the converging fiow. The dominant contribution is in oq (negative) which varies linearly 
with the rotation rate - figure 9(d) - while the axisymmetry imposes (ai) = (03) = 0. 
When one probes Bz, the dominant contribution comes from the differential rotation whose 
effect is not felt on the rotation axis, so that (ao) = 0. The dominant contribution is in 
(ai), which varies linearly with Q (figure 9(b)). 

• Fluctuations of the polynomial coefficients 

Let us first discuss the case of the fiuctuations in the Omega effect, i.e. for measurements 
of the toroidal field induced by differential rotation when the applied field is axial. The 
time variations of the leading coefficient ai are shown in figure 10(a). One observes very 
strong fiuctuations, although the probability density function is quasi- Gaussian - Fig. 10(b). 
One finds ai^i.jns/('^i) = 114% at a (counter) rotation rate f2 = 10 Hz. The spectrum in 
figure 10(c) indicates that the fiuctuations in the profile have a long-time behavior (the 
slope of the spectrum in the low frequency range is close to -0.7). One thus observes that 
the fiuctuations in the induction profile correspond to a slow process, compared to the 
turbulence fast small scales and also compared to the disk turn-over time. The fact that the 
entire induction profile changes in time is clear in figures 10(d,e) which show the correlation 
between the coefficients Oq and ai. One observes that Oq and ai are anti-correlated, with 
ao ~ only for small deviations of ai about its mean. In most configurations Oq is non-zero. 
As at small Rm the induction mirrors the evolution of the velocity, we conclude that the 
fiow has a slow dynamics with strong deviations from the mean von Karman geometry. 

We now compare the fiuctuations in time for the leading polynomial coefficient and each 
choice of (i?o,i,-Bj) - figure 11. Two kinds of behavior are evidenced. When {i,j) = {z,x) 
or {i,j) = {x,z), i.e. as one probes the toroidal fiow (infiuence of differential rotation) 
the fluctuations are large. When the induction probes the poloidal component (stagnation 
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FIG. 10: Fluctuations of the polynomial coefficients, for the Omega effect: a magnetic field i^o.z = 
24 G is applied along the rotation axis, and the induction profile Bg = is measured. Counter- 
rotation of the disks at = 10 Hz. (a) time evolution of the polynamial coefficient of order 1; 
(b) corresponding centered probability density function; (c) time spectrum; (d) and (e): cross- 
correlation and joint probability density function for the coefficients ao{t) and ai(t). 

point), i.e. for = {z,z) or = {x,x), the fluctuations are reduced by a factor 3. 
Note that in figure 10 the coefficients are compared in gauss, so that the ai values have been 
multiplied by a length chosen as half the cylinder radius. This choice is justified because in 
the profiles the induction is maximum at R/2 in the corresponding case. In addition, R/2 
is also close to the diffusive length. Another feature is the change in the ratio of the rms 
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fluctuation amplitude to the mean. In the case of an applied field parallel to the rotation 
axis, we had ai^rms/(ai) = 114% at a (counter) rotation rate f2 = 10 Hz for induction 
resulting from the differential rotation. For the induction due to the pumping motion 
towards each disks, one has ao,rms/ (c^o) = 20%. Note that a 20% fluctuation level is what 
one observes for induction effects in the case of a single rotating disk One thus finds 
that fluctuations in the induction are larger when associated to the toroidal flow compared 
to the poloidal velocity component. This is confirmed in figure ll(c,d) where the applied 
field is transverse to the rotation axis. There Oq for the compression effect fluctuates with 
a level ao,rms/(a.o) = 50%, while ai for the BC-effect has ai,rms/(a.i) = 160%. In all cases, 
the PDFs of the time fluctuations of the coefficients are quasi-Gaussian. 



(a) Axial field : £2 effect 
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FIG. 11: Fluctuations of the polynomial coefficients, for four orientations of the applied and 
measured components, {Bo.i,Bj). (a) = {z,x); (b) = {z,z); (c) = {x,z); (d) 

= (x,x). Counter-rotation at $7 = 10 Hz. The solid red line corresponds to the mean value 
and the dashed line to a level of fluctuation equal to 20% of the mean. 
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• Very slow quasi-periodic modes 

We must mention that in figure 11(d), the profiles were high-pass filtered at a frequency 
above 3 Hz (for a disk rotation rate of 10 Hz) before processing . The reason is that there is 
a quasi-periodic evolution in the induction profile, as can be observed in the signal displayed 
in figure 12 (and previously in Figure 4(b)). The unfiltered time variations of ao(^) are 
shown in figure 12(a), and its corresponding time spectrum and probability density function 
in figures 12(b,c). The spectrum has a marked peak at a frequency / = 0.1 Hz ~ ^2/100, with 
the quasi-periodic evolution also reflected in the bimodal shape of the probability density 
function. 
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FIG. 12: Transverse field and compression effect : applied field Bq^x and induced component B^. 
(a) time variation of tfie leading polynomial coefficient oq (low-pass filtered at frequencies lower 
than 3 Hz); (b, c) time spectrum and probability density function of ao{t). Measurement with 
counter-rotating disks at 17 = 10 Hz. 



In this configuration, the external magnetic field is transverse. The induced component 
probed is parallel to the applied field, and results from the compression of the applied field 
lines by the poloidal flow converging in the median plane. This effect strongly depends 
on the exact location of the shear layer Its oscillations with respect to the mid-plane 
would generate oscillations as in Fig. 12(a). We note that these slow oscillations are 
reminiscent of the global instabilities of the median shear layer, previously discovered and 
studied in water prototypes by CEA-Saclay group jl^. |27|. 
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IV. DISCUSSION AND CONCLUDING REMARKS 



In the last case discussed above, the fluctuations are attributed to the existence of large 
scale (vortical) structures, possibly due to the roll-up of the shear layer in the von Karman 
flow. When these structures exists, it is obvious that the flow is no longer in the averaged 
configuration pictured in figure 1(a). We would like to show now, that even if one discards 
the influence of these large scale coherent structures (as we have done by high-pass filtering 
the data), the inherent turbulence of high Reynolds number von Karman flows is such 
that the fluctuations in the instantaneous flow geometry are very large. The instantaneous 
flow configurations differ from the mean flow, not only in regards to the amplitude of the 
toroidal and poloidal velocity components, but also in regards to the symmetries and to the 
flow overall structure. 

We first comment again about the fluctuations, and we study in greater details the 
induction profile for the Omega effect. In this case the induction is less sensitive to 
displacements of the mid-plane shear layer because the differential rotation is found to 
be rather uniform in the center of the flow For {Bqz,Bj.), the coefficient should 

vanish because of axisymmetry. Its actual fluctuations in time are shown in figure 13. The 
mean value (ao) is equal to 0.5 gauss; although one expects (oq) = 0, the recorded mean 
value is really within the precision of our polynomial analysis (limited by the number of 
elements in the array probe). On the other hand, the rms amplitude of fluctuations for 
ao is ao,rms = 1-25 G (for an applied field equal to 24 G), well above its mean. This is a 
significant variation, even when compared to the mean induction {{ai)R/2 = —2.5 gauss). 

Let us now discuss changes in the flow geometry. We compute how long the flow remains 
in the neighbourhood of its mean configuration, by imposing that both ai and Oq take values 
within a standard deviation of their mean. In this case, one finds again that no longer than 
50% of the time is spend in the mean von Karman geometry. We stress that: (i) a similar 
conclusion is reached if one does the estimation from induction effects generated by the 
presence of the stagnation point (stretching of the axially applied field lines or compression 
of a transverse field); (ii) the large changes in the polynomial coefficients are related to the 
long-time dynamics of the flow {i.e. to frequencies lower than that of energy injection into 
the flow). 
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FIG. 13: Q effect, for a counter rotation at 10 Hz. Time fluctuations of the oq coefficient (shifted 
by 0.5 G (see text) so that (oq) = in the plot). 

At the low magnetic Reynolds numbers probed in this Gallium experiment, the above 
effects can be attributed to the hydrodynamics of the flow. Note that such is not neces- 
sarily the case for sodium experiments in which non-linear induction processes Q, 13] or 
dynamo 0, 0] take place. The understanding of such long time dynamics in confined tur- 
bulent fiows is a challenge. It has now been reported in several experiments (e.g. recently 
in Rayleigh-Benard convection [2^) but its understanding is still elusive. For instance no 
known argument gives the time scale of the slowest motion. 

These slow global changes in the geometry of the fiow may not favor dynamo action. 
First because all configurations may not be consistent with self generation - in the sense 
that the velocity field U(r, t fixed) may not lead to a positive growth rate when inserted 
in the induction equation. Second, in order for the the magnetic field to grow, the fiow 
must be maintained for many kinetic advection times; indeed one has for the magnetic 
diffusion time tm = R^/\ ~ ^mTNL, with tnl = R/U the time scale of the forcing. 
Altogether these arguments indicate that a stable flow conflguration is desirable for the 
self-generation of a stationary dynamo. For instance, for flows generated by the rotation 
of only one disk with the other kept at rest, we have observed that the fluctuations are 
much less — cf. Figure 14. The flow spends most of its time in the Siti conflguration 
imposed on average by the driving. As a result, depending on the forcing and large 
scale hydrodynamic evolution of the flow, one may have to be cautious with approaches 
that estimate the dynamo threshold from mean flow geometries. This procedure may be 
valid for purely helical flows in which we have not detected strong fluctuations about 
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the mean; it certainly did yield a correct estimate of the onset of the Riga dynamo j5[. 
However, it may not be the case for the flow generated by counter-rotation in the von 
Karman geometry. The observed slow dynamics is associated with important changes 
in the flow topology, and a mean field kinematic simulation may underestimate the threshold. 








0.5 1 1.5 2 2.5 3 3.5 



E(G) 

FIG. 14: Probability density function of the fluctuations of the £.2 distance Eif) to the mean 
induction profile. (a,b) flow generated by counter rotating disks at 17 = 10 Hz, axial applied fleld, 
(a) : induction due to differential rotation, (b) due to stretching; (c) flow generated by the rotation 
of one disk only at 17 = 10 Hz, transverse applied field. 
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